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ABSTRACT 

We make predictions for the metallicity of diffuse stellar components in systems 
ranging from smal l spiral galaxies to rich galaxy clusters. We extend the formalism of 
iPurcell et al.l (|2007i) . in which diffuse stellar mass is produced via galaxy disruption, 
and we convolve this result with the observed mass- metallicity relation for galaxies in 
order to analyze the chemical abundance of intrahalo light (IHL) in host halos with 
virial mass lO^^'^M© < Mhost < lO^^M©. We predict a steep rise of roughly two dex 
in IHL metallicity from the scales of small to large spiral galaxies. In terms of the total 
dynamical mass Mhost of the host systems under consideration, we predict diffuse light 
metallicities ranging from Zjhl< — 2.5 for Mhost ^ IO^^Mq, to Zjhl ^ —1.0 for 
Mhost ~ IO-'^^Mq. In larger systems, we predict a more shallow rise in this trend with 
ZiHL - -0.4 for Mhost - IOI^Mq, increasing to Zjhl - 0-1 for Mhost ~ lO^^Mg. 
This behavior is coincident with a narrowing of the intrahalo metallicity distribution 
as host mass increases. The observable distinction in surface brightness between old, 
metal-poor IHL stars and more metal-rich, dynamically-younger tidal streams is of 
crucial importance when estimating the chemical abundance of an intrahalo population 
with multiple origins. 
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1 INTRODUCTION 

Diffuse stellar components are ubiquitous phenomena 
that include the intracluster light at the centers of 
rich galaxy clusters as well as the diffuse stellar halos 
that surround individual galaxies like the Milky Way. 
This diffuse luminosity comprises < 1 — 5% of t he to- 
tal luminosity of spiral galaxies (fMorrison et al.l 
Chiba fc Beers' '2000*; 'Yanny et al.' '2000'; 'ivezic et al 
Sieg cl ct al. 2002; Irwin ct al. 2005; Guliathaku rta et al.l 
20051 : IKalirai et al.l 120061 : IChapman et al.l bOOfiP and as 
much as ~ 10 — 40% of the luminosi ty of rich galax; 
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clusters lICalcaneo-Roldan et al.' 2000; Lin fc 
Feldmeier et al.ll2004l; 'Mihos et al. 2005i; iZibetti et al.l 
Krick et al. 2006; Scigar ct al. 200,6). Whether in clusters 
or surrounding individual galaxies, we refer to this dif- 
fuse material as "intra-halo light" (IHL). In the prevail- 
ing theory of structure formation, objects are built hier- 
archically through sequential mergers of smaller objects. 
This provides a natural means of IHL p roduction via the 
disruption of merging galaxies ( IGallaghcr fc Ostrikei] 1972 : 
Merritt"l983':'Bvrd fc Vahonen' 1990; Dubinski et avhoOS; 
Gnedin 2003; Mihos 2004; Murantc ct al. 2004; Lin fc Mohr 
20041 ') . This mechanism includes the emergence of diffuse 



stellar halos about galaxies such as the Milky Way, resulting 
from the merger and dest r uction of dim, dwarf-s i zed galax- 
ies jSearle fc ZinnI Il978l; Ijohnston et all [l996l: iJohnstonI 
199a lHelmi et al.lll999l:lBullock et al.ll200ll:[johnston et all 



20051 : iRobertson et all |2005| : 



20011: IBuUock fc Johnston 

De Lucia fc Helmil bOO^T making the IHL a potentially rich 
testbed for structure and galaxy formation theories. In a 
previous paper, we made predictions for the relative promi- 
nence of the IHL as a function of both the total mass of the 
system (including luminous and dark matter) and the lu- 
minosity of the brightest object in each syste m for systems 
rangi ng from dwarf galaxies to rich clusters l|Purcell et al.l 
l2007h . In this paper we seek to add to the discriminatory 
power of the IHL as a test of galaxy formation models by 
making predictions for the metallicities of IHL components 
in such systems. 

Observations of intracluster light in Abell 3888 indicate 
that the stellar population composing it has a metallicity 
distribution function peaking near Z ~ 1.0 in the outer 
regions of the di ffuse component, and 0.2 < Z < 0.5 in 
the inner region (iKrick et al.ll200a ). although there is large 
uncertainty in both values stemming from chemical evo- 
lution systematics and color analysis. Conservatively, we 
can say that the chemical abundance of intrahalo stellar 
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Figure 1. The global mass-to-light M/Lc ratio as a function 
of host halo mass Mhaio- The so l id lin e shows the M/Lc result 
obtained by Ivan den Bosch et alj ll2007t) for the WMAPl cosmo- 
logical parameters. For comparison with the model of PBZ07, the 
dotted line repres ents that work's fiducial M/Lc as drawn from 
lYang et alJ l|2003h . 



material in cluster-sized host systems is at least one or- 
der of magnitude greater than that of the diffuse stellar 
components around large spirals. M31's pressure-supported 
halo has an approxima te iron a bundance of [Fe/H] ~ — 1.4 
jChapman et all |20od see also iKalirai et alJbOOSl ). and it 
resembles the Milky Way halo both kinematically and struc- 
turally. The Milky Way's halo stars are also metal- poor at 
around [Fe/H] ~ -1.5 {e.g., iRvan fc Norr3 119911 1. How- 
ever, M33's halo has proven elusive as there is only ten- 
tative evidence for a non-rotating, power-law diff use com- 
ponent at roughly [Fe/H ] ~ —1.5 (|Fergusonll2007l . see also 
iMcConnachie et al.ll200a for spectroscopic field analysis). 

Numerical simulations indicate that the shredding of 
infalling galaxies during the hierarchical merging process 
of host halo form ation is a prevalent contributor to diffuse 
light in c lusters (^Willman et~al] l2004l : iRudick et all l2006l : 
ISommer^L arscn 2006). There is often sufficient stellar mass 
in this component to dominate th e light emitted by the 
bright central galaxy in the cluster |Conrov et aLll2007l . see 
iGonzalez et al.li2005i : ISeiear et al.ll20oi for observational ev- 
idence supporting this prediction). Likewise, simulations of 
Galaxy formation support the accumulation of a spheroidal 
stellar halo via the accretion and disruption of dwarf-type 
satellite galaxies during the construction of the primary 
system toiemand et al.ll2005l : iRead et al ] |2006l : iFont et al.l 
l2006l : lAbadi et al.ll2006l i. 

Our aim in this paper is to predict the metallicity of 
intrahalo hght as a function of the total mass of the system. 
We explore the behavior of any trends in metallicity as a 
function of properties of disrupted progenitor systems, in or- 
der to discriminate bright diffuse streams from dynamically- 



older and more well-mixed tidal debris. This debris is created 
from an enormous dynamic range in halo mass. In addition, 
the system-to-system scatter in the fraction of luminosity in 
IHL and IHL metallicities is significant, indicating that large 
samples are needed in order to study IHL metallicity trends 
as a function of mass. These factors make direct simulation 
of IHL production a computationally costly endeavor that 
is impractical for our purposes, particularly given the inher- 
ent uncertainties in numerical treatments of baryonic pro- 
cesses. As an alternative, we rely on the analytic treatment 
prescribed by IZentner et al.l (|2005l . see below) in order to 
follow the coalescence of small systems into larger compos- 
ite systems. This approach models the merging rates and 
subsequent interactions, including heating and disruption, 
of systems in an approximate way but with the benefit of 
being a computationally inexpensive approach that can be 
used to minimize statistical (not systematic) uncertainties 
in model predictions. We assume that the stellar material in 
a merging galaxy will be diffused into the intrahalo light of 
the larger system when the galaxy's parent subhalo is signif- 
icantly affected by tidal and heating processes. Direct treat- 
ments of baryonic processes and star formation are uncer- 
tain at present. As such, we employ a "bootstrap" approach 
that feeds the observed mass-metallicity relation for galaxies 
into our calculation. After determining the stellar contents 
of disrupted, mergin g halos according t o the empirically- 
constrained model of lPurcell et al.l l|2007l . see below), we as- 
sign these stars a metallicity based on the observed mass- 
metallicity relation. The metallicities of all stars distributed 
throughout the diffuse light by all merging and disrupted 
systems are used to obtain the oxygen abundance [0/-ff]iHL 
of diffuse intrahalo light as a function of the total mass of 
the system Mhost- At this time, detailed and robust predic- 
tions remain beyond the reach of theoretical modelling, so 
we em phasize (as we did in our previous study. IPurcell et al.l 
l2007f ) that our goal is to highlight gross trends that should 
be expected in hierarchical models on very general grounds. 

In the next section, we outline our model for IHL metal- 
licities. In §[3l we present our results and discuss their impli- 
cations in the context of observational benchmarks in chem- 
ical abundance measurements of galactic stellar halos as 
well as intragroup/intracluster luminosity. Throughout this 
work, we adopt a ACDM cosmological model with h — 0.7, 
i^m = 1 — = 0.3, and a primordial power spectrum which 
is scale-invariant, n = 1, and normalized to erg = 0.9. 



2 METHODS 

We investigate the metallicity of tidally-disrupted, diffuse 
stellar material as a function of host dark matter halo mass. 
In our terminology, we refer to the largest dark matter halo 
surrounding a group or cluster of galaxies as the "host" dark 
matter halo and the smaller halos that are themselves con- 
tained within the host but self-bound syste ms as "subha- 
los." Our study expands upon the model of iPurcell et al.l 
(|2007l . hereafter PBZ07), which makes statistical predictions 
for the stellar mass fraction expected in a diffuse, luminous 
intrahalo component, over a wide range of host halo mass 
(lO" '^ Mq < Mhost < 10^^ Mq). We give a brief outline of 
the relevant features of this model here. A more complete 
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discussion of our modeling techniques and definitions can be 
found in PBZ07. 

Out treatme nt utili zes the halo formation model 
of Zentner et al. (120051 . for an earlier version see 
IZentner fc BuUocM l2003h in order to track the accretion, 
evolution, and possible dest ruction of dark matt er subha- 
los analytically. The model of lZentner et al.l l|2005h produces 
mass accretion histories over a wide range in host halo mass 
acco rding to the extended Press-Schechte r (EPS) formal- 
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ism llBond et al.lll99ll: iLacev fc ColdflQQ^ . for a recent re- 
view see IZentneil'2007h^ using th e spec ific implementation 
described bv Somcrvill e fc KolattI l| 19991 ). and then uses an- 
alytic prescriptions for dynamical friction and mass loss due 
to interactions within the larger host potential to evolve the 
population of accreted subhalos. This model is very well- 
tested against dissipationless cosmological simulations, re- 
producing a variety of subhalo statistics predicted by direct 
numeri cal techniques accurat ely. We encourage the reader to 
refer to lZentner et al.l (|2005l ) for more detail on this aspect 
of our model. 

Applying observational constraints from galaxy survey 
data, PBZ07 assign subhalo luminosit ies according to the 
results obtained by lYang et al.l (|2003l . see below) in their 
analysis of the conditional luminosity function of galax- 
ies in the Two-Degree Field Galaxy Redshift Survey (2dF- 
GRS). The model of PBZ07 then implements a critical 
mass-loss threshold at which orbiting subhalos are consid- 
ered to be completely tidally shredded when their circu- 
lar velocities fall below a critical fraction of their maxi- 
mum circular velocity they had upon accretion into the host, 
/crit = Vcrit/Vmax{tacc)- A sccoud free parameter a describes 
the star formation history of the subhalo, in order to deter- 
mine its stellar-mass content upon accretion and the subse- 
quent truncation of star formation via ram pressure strip- 
ping induced by the host halo's hot gas. Both parameters 
are tuned such that contemporary populations of surviving 
galaxies are recovered by the fiducial model across the host 
mass scales of interest, and it is instructive to note that the 
results are insensitive to dramatic changes in these values, 
which are set to fait =0.6 and a = 1. We ignore contribu- 
tions made to a system's diffuse light by surviving satellites, 
as simulations indicate that subhalos are typically disrupted 
shortly after t heir stars become tidally st ripped in any sub- 
stantial sense (|Bullock fc Johnston|[2005h . 

The global variance in the mass-to-light r atio of dark 
halos as a function of host mass, as obtained bv lYang et al.l 
and discussed in § 12.11 convolved with the EPS for- 
malism's prediction that host halos are chiefly built by the 
accretion of subhalos roughly a tenth as massive, yield a 
clear correlation between the relative amount of intrahalo 
light in a system and the total mass of the system. This pre- 
dicted correlation spans three orders of magnitude in mass 
from galactic scales (/ihl ~ 1 — 5% for Muost ~ 10^^ Mq) 
to those of galaxy clusters (/ihl ^ 10 — 40% for Almos t ~ 
lO^'^Mo), and we refer the reader to lPurcell et al.l l|2007h for 
a complete description of the method and prior results. 



2.1 Intrahalo Light Model Update 

An ingredient leading to the correlation between IHL frac- 
tion and host mass is the aforementioned variation of mass- 
to-light ratio [M/Lc) with halo mass, which is necessary 
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Figure 2. The predicted IHL fraction as a function of host halo 
virial mass ranging from W^'^-^Mq < Muost < W^^Mq. The dia- 
monds denote the mean of the distribution at fixed mass based on 
1000 realizations of our analytic model, and the solid line marks 
the median. The light shaded region encapsulates the 95% range 
of the distribution at fixed mass, and the dark shaded region 
contains 68% of the distribution. The upper axes show the cen- 
tral galaxy (BHG) luminosity as derive d from the M/^aat/ Lbhg 
mapping of lvan den Bosch et al. I ||2007| '). 



in the prevailing structure formation model in order to 
reproduce empirical measurements of the g alaxy luminos- 



ity function and clustering s t atistics (e.g., White fc ReesI 



1978i; iKauffmann et all Il993l: ISomerviUe fc Primackl Il999l : 



Tinker et al.1 l2005l : ICoorav fc Milosavlieyid I20q5^ . PBZ07 



adopt the particular result of lYang et al. n|2003h . in which 



the authors constrain the conditional luminosity function 
(CLE) of galaxies in the Two-Degree Field Galaxy Redshift 
Survey (2dFGRS) and infer the characteristic B-band lumi- 
nosity, Lc(M), for the brightest central galaxy sitting in a 
host h alo of virial mass M. In a n ew analysis by the same 
group, Ivan den Bosch et al.1 (|20o3) investigate variation in 
their CLE results as a function of the cosmological parame- 
ters measured by WMAP0, and introduce an updated theo- 
retical model, which is integrated over the light-cone interval 
defined by the upper and lower redshift limits given by 2dF- 
GRS data, and accounts for the scale-dependence of halo 
bias. 

In Figure [J we sh o w the mass-to-light ratio derived 



bv lvan den Bosch et al. (^2007^ and comp are it to the pre- 
vious CLE analysis of Yang et "al] (|2003l ). which was used 



in the fiducia l model of PBZ07. I n this paper, we update 
the model of iPurceU et al.1 ||2007|) to refiect the more re- 
ce nt analysis of van den Bosch et al.1 l|2007h . The results 
of lvan den Bosch et al.1 (|2007l i indicate that in a WMAPl 
cosmology galaxy formation occurs most efficiently in ha- 
los of host mass M ~ 10"-^ M©, where mass-to-light ra- 
tio is a minimum. Otherwise, M/Lc increases as a power- 



^ In this paper, we will refer to three cosmologies corresponding 
to WMAP data releases: WMAPl with = 0.3 and as = 0.9, 
WMAP3 with n„ = 0.23 and erg = 0.74, and WMAP5 with 
n„ = 0.258 and as = 0.796. 
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Figure 3. The fiducial, redshift-dependent mass-metallicity re- 
lation used to assign oxygen abundances to disrupted satellite 
galaxies produced by the model of PBZ07, in order to assess the 
ch emical content of int rahalo luminosity. Drawn from the model 
of lSavaglio et al.l ||2005| ') and formulated explicitly in Equation 1, 
the maroon lines represent case 1 of the M* — Z relation at z =0, 
1, and 2 (solid lines of high, medium, and iow thicknesses, respec- 
tively), and case 2 is plotted in orange dashed lines at the same 
redshifts, representing the original relation shifted by 0.3 dex to- 
wards lower metallicity, to reflect selection effects discussed in 
i|2.2l For comparison, we overplot with black squares data from 
the m ass-metallicity relation at z ^ 2, as derived by l&b et al.l 



law with M/Lc oc where 7 — —2.32 for systems with 
M < 10" '^ ft-^Mo, and 7 = 0.27 for M>10" '^ fe-^Mp). By 
compa rison, the M/Lc ratio derived bv lvan den Bosch et al.l 
l|2007h for the WMAP3 cosmology is 25-45% lower across 
the spectrum of host mass, with the most pronounced dif- 
ference occurring at the minimum near Mhost ~ lO^^'^M©. 
We adopt the WMAPl M/Lc mapping in order to remain 
consistent with the cosmological parameters used in our pre- 
vious work, but IHL fractions produced by our model in a 
WMAP3 cosmology differ by less than 1% from the WMAPl 
result on all mass scales, and we note that the intermediate 
WMAP5 parameters therefore also yield indistinguishable 
results for the IHL fraction. In addition, we draw our lumi- 
nosity values according to Gaussian distributions of M/Lc 
at fixed host mass, with a constant a — 0.4 dex variance that 
roughly approxima tes the scatter in the proba bility distri- 
butions obtained bv lvan den Bosch et al.l (|2007| ). though our 
results are relatively robust to this choice. 

A less significant alteration to the PBZ07 model in- 
volves additional recursion levels applied to the EPS merger- 
tree algorithm, so that we now model the orbital evolution 
of subhalos-of-subhalos at each host mass, allowing us to 
probe diffuse light already present in an infalling system. 
Numerical simulations suggest this pre-processed IHL may 
be a dominant contributor to the mild increase in the diffuse 
light fr action as a function o f incre a sing host mass on clu ster 
scales (|Murante et al. 1 12004 l2007l : [Monaco et al.l [2006! ) . as 
observed by surveys testing the evolution of intracluster light 



with cluster richness {e.g., IZibetti et ahllioosl ). In Figure [2j 
we present updated results for the relative fraction of lumi- 
nosity in diffuse IHL compared to the total luminosity in the 
system, which includes the contributions from the brightest 
halo galaxy (BHG) in the center as well as destroyed (IHL) 
and surviving satellites (Sats), as a function of host halo 
mass. We predict a positive correlation between halo mass 
and IHL fraction spanning more than an order of magnitude 
in IHL fraction, with the dominant factor in the variation be- 
tween this model and the result of PBZ07 being the updated 
slopes and mini mum value in the beh a vior o f M / Lc[Mhost 
as presented by Ivan den Bosch et al.l (|2007l l. This predic- 
tion yields distributions of intrahalo light roughly consis- 
tent with galactic stellar halo luminosities, /ihl ^ 1 — 5% 
for Mhost ^ 10^^ M0, as well as intracluster observations of 
diffuse fight, /ihl ~ 10 - 40% for Mhos 



2.2 Assigning Metallicities 

Having predicted tfie relative amount of intrahalo stellar 
material in systems ranging in mass from dwarf galaxies 
up to rich galaxy clusters, we can also determine the ex- 
pected chemical abundances in these diffuse components 
by considering the well-known correlation between luminos- 
ity a nd metallicity (L-Z relation) in star- forming ga laxies 
[e.g., iGarnett fc Shield J 19871 : iBrodie fc Huchralll99ll ). The 
L-Z relati on, which holds over ten magnitudes in optical lu- 
minosity (jZaritskv et al.l [19941 ). is generally thought to be 
the result of an underlying correlation between stellar mass 
and metallicity for galaxies. Tfiese two properties are fun- 
damental relics of galaxy formation and evolution, because 
they refiect the amount of gas tied up in a system's stars 
and the efficiency of gas recycling and retention processes in 
that system, respectively. 

Resolving the age-metallicity-reddening degeneracy for 
a particular stellar population translates a galaxy's luminos- 
ity into stellar mass, M*. The — Z correlation is osten- 
sibly a consequence of the interplay between a host halo's 
total baryonic mass and the effective yield of the system; 
larger galaxies being less susceptible (having higher stellar 
mass quotients) to outflow processes that decrease the quan- 
tities of m etals present for sub sequent generations of star 
formation ijTremonti et al.ll20o3 ). 

Extensive surveys of star-forming galaxies have con- 
strained the mass-metallicity relation at various redshifts. 
Using oxygen gas-phase abundances derived from optical 
nebular emission lines in gala xy data drawn fr o m the Sloan 
Digital Sky Survey (SDSS), iTremonti et"aLl (|2004l ) infer 
an Mi, — Z relationship for the local universe (2<0.1), 
which is analytically well-fitted by a second-order polyno- 
mial in log(M«). This trend applies for galaxies with stel- 
lar massesbetween 10* M0 < Mi, < lO^ '^M©, and the work 
of I Lee et al] l|2006l ) extends the correlation into the low- 
mass regime {IQ^Mq < Mi, < lO^'^M©) by measuring 4.5 fim 
luminosities for nearby dwarf irregular galaxies observed 
by near-infrared instrumentation aboard the Spitzer Space 
Telescope. For the mass range in which the two determi- 
nations overlap, the SDSS relation is roughly 0.3 dex more 
metal-rich than that derived for the dwarf-galaxy data, an 
offset that may be partially explained by two selection ef- 
fects: SDSS fiber spectra preferentially sample the innermost 
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Figure 4. The oxygen abundance, relative to the solar value, for diffuse intrahalo stars in systems of mass IO^'^'^Mq < Mhost < 10^^ Mq, 
in both cases of the model. The doubly thick solid line in each panel represents the mean of the [0/-ff]iHL distribution at each host mass. 
The thick solid and thin solid lines denote the boundaries containing 68% and 95% of th e distribu t ion, r espectively. The observational 
estimate for the Milky Wa y stellar halo is shown with a diamond, Z ~ — 1.5 as reported bv lFergusonI ||2007|). w i th an added o-enrichment 
of +0.4 jMcWiUiamlllQQTt) . The Galactic virial mass is taken to be 2.43 X lO^^M© as derived by I Li fc Whitj ll200d) . 



regions of spiral galaxies w hich may in actualit y have radial 
metallicity gradients {e.g., IZaritskv et al.|[l994l ). and oxygen 
abundances derived for SDSS data use bright-line estimation 
methods (in the absence of the temperature-sensitive [O III] 
emission line) which overestimate the true metallicity val- 
ues b y '^0.3 dex at roughly solar abundance a nd greater 
{e.g., iKennicutt et ahlbood : Ifiresolin et al.ll2004h . 

There is evidence that the mass-metallicity relation has 
evolved with redshift. Specifically, galaxies with a given stel- 
lar mass at high redshift tend to have lower oxygen abun- 
dances than_their_Jocal_coii^^ same stellar 
mass (|Savagho et~aLll2005l : ferb et al. ' '20061. Armed with a 
mass accretion history for disrupted subhalos, we can assign 
their stellar mass values using the adapted PBZ07 model 
and determine the appropriate oxygen abundance of the re- 
sultant IHL component by applying an AI-^ — Z relation that 
evolves with redshift to match each satellite's epoch of de- 
st ruction. To accomplis h this goal we adopt the prescription 
of lSavaglio et al.l (|2005l . hereafter 805), in which the authors 
investigate stellar mass and metallicity in galaxies at redshift 
0.4 < 2 < 1.0 from the Gemini Deep Deep Survey (GDDS) 
and Canada- France Redshift Survey (CFRS), in order to de- 
rive the form by which the Mi, — Z relation changes as a func- 
tion of cos mic time. This model t ranslates the M* — Z poly- 
nomial of iTremonti eraP (j2004h into one consistent with 
their choice of stellar IMF and metallicity calibration, and 
extends the relation in redshift. It is worth noting here that 
we determine stellar metallicity via the time-dependent for- 
mulation describing oxygen gas-phase abundance, such that 
the chemical composition of our stars will reflect the metal- 



licity of the gas from which they were formed, instead of 
the systematically richer gas-phase metallicity that would 
be measured at the present day. 

To illustrate the importance of the selection-effect issues 
discussed above, we explore two mass-metallicity relations in 
our work. We present predictions for IHL metallicity based 
on each of these mass-metallicity relations. The two M* — Z 
relations represent relatively high- and low-metallicities re- 
spectively and we expect that the predictions from these 
two relations bracket the true value. In case 1, we adopt the 
model of SOS, and in case 2 we decrease the oxygen abun- 
dance by 0.3 dex across all stell ar masses as mo tivated by 
the offset found in the analysis of lLee et al.l (|2006t ). For each 
destroyed satellite galaxy in our analytic formalism, we draw 
from a Gaussian distribution of metallicities centered on the 
mean M* — Z relation and with a dispersion constant over 
host halo mass and equal to cr = 0.2 dex as derived by SOS 
for their sample of galaxies si z ^ 0.7. In Figure O we plot 
the fiducial mean mass-metallicity relation at 2 = 0, 1, and 
2 for both cases. Throughout this work we adopt the stan- 
dard units for oxygen metallicity in g alaxies, 12-|-log(0/H ) 
in which the solar abundance is 8.66 (jAsplund et al.ll2004l l. 
The functional form of our case 1 as drawn from SOS is, for 
a given stellar mass M-t, and Hubble time tu'. 

12 + log(0/H) = - 7.S9 + 2.S31S log Af* - 0.09649 log^M* 
-)- S.1733 log tu - 0.3944 log^tn 
- 0.403 log tH log M* (1) 

For this prescription, in order to reflect the buildup of both 
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Figure 5. The [0/-ff]iHL distributions as a function of host mass for diffuse stellar populations above or below a threshold in surface 
brightness. In the left panel, we select dynamically older intrahalo light by only showing the IHL metallicity for disrupted subhalos having 
dimmed to fig > 28 mag/arcsec^, while the center panel includes diffuse material less bright than a slightly higher threshold of /ig > 27 
mag/arcsec'^. By contrast, the right panel selects only dynamically young IHL features that remain brighter than fig = 27 mag/arcsec'^. 
The thick solid lines denote the mean of the distribution at fixed host mass, while the thin solid lines are contours bracketing 68% of 
each distribution. 



stellar mass and the associated chemical composition of sub- 
halo stars, we adopt the tg at which the satellite galaxy has 
formed half of the stellar mass it has at the time of accretion. 
I.e., the appropriate A'U{tg) = 0.5M*(tacc). As in PBZ07, 
we have assumed for simplicity that star formation is trun- 
cated upon the subhalo's accretion onto the host halo, via 
processes of gas dynamics such as ram pressure stripping or 
strangulation. 



Remnants of relatively recent accretion events will nat- 
urally be more metal-enriched than the diffuse Ught from an- 
cient merger remnants, which give rise to a more spatially- 
uniform, nearly-virialized component of diffuse luminosity 
that should be more metal-poor than the tidal streams 
that are the cohe rent remnants of recent accretion events 
ijPont et al.l l2006l ) . Our formalism defines subhalo disrup- 
tion by mass loss alone and we do not trace the subsequent 
evolution of the liberated material in a self-consistent man- 
ner. Consequently, we must introduce additional modeling 
in order to distinguish bright tidal features from older relics 
of galaxy disruption. In their paper on the interpretation of 
diffuse streams as the result of satellite disruption around 
galaxies in a ACDM universe. Ijohnston et al.| (|200ll 'l derive 
the form for the time evolution of a debris trail's surface 
brightness, as a function of the progenitor satellite's mass, 
the host halo mass, and the galactocentric radius of the 
satellite's orbit. Their analysis is based on the destruction 
of satellites on circular orbits, but we utilize this approxi- 
mate form in order to test the viability of separating the IHL 
metallicity along thresholds in dynamical age, which could 
be useful in studies of Local Group diffuse features such as 
the giant southern stream of M31, as we discuss in ^ 



3 RESULTS 

We make predictions for the oxygen metallicities of diffuse 
intrahalo stellar mass embedded in host halos with virial 
mass IO^^'^Mq < Mhost < W^^Mq, as shown explicitly in 
Figure m In case 1, corresponding to the redshift- dependent 
Mi, — Z relation obtained by SOS, our analysis predicts a 
mean [0/-ff]iHL of ~ —0.2 — 0.2 with respect to the so- 
lar abundance, for diffuse intracluster luminosity present 
in halos with mass larger than Mhost ~ IO^^Mq, while 
the case 2 metallicity is more sub-solar at those scales, 
[0/H]iHL ~ —0.5 — 0.0 as expected due to the differ- 
ence between the normalizations of the two fiducial models. 
On galactic scales, we expect the abundance of the stel- 
lar halo to rise quickly as a function of host mass, increas- 
ing from [0/H]mL ~ —2 to [0/H]iul ~ — 1 in the range 
1O"'^M0 < Mhost < IO^^Mq. IHL metallicities in both cases 
evolve much more slowly in the regime of intragroup and 
intracluster light, the trend growing more shallow for host 
halos larger than Mhost ^ lO^'^'^ and rising weakly from 
[0/H]iuL ~ —0.5 to roughly solar metallicities, inducing a 
more distinct separation between case 1 and case 2, whose 
[0/H]mh distributions bracket the solar abundance on the 
largest scales investigated. 

The distribution of IHL metallicity at fixed host mass 
is quite wide for stellar halos, with variance az ~ 0.5 dex 
at Galactic-mass scales, and is significantly more narrow 
for intragroup and intracluster light, where az ~ 0.1 dex 
and the average chemical abundance of diffuse stellar mass 
reaches a plateau. The normalization of this plateau varies 
non-trivially as we consider imposing thresholds in surface 
brightness, corresponding to the selection by dynamical age 
of destroyed satellite galaxies in order to distinguish the old. 
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Figure 6. The difference between tfie metallicity of a liost halo's 
central galaxy, Zji^at^ S'lid the metallicity of the intrahalo compo- 
nent, ZjHL- For each case, the thick solid lines denote the mean 
of the distribution at fixed host mass, while the thin solid lines 
are contours bracketing 68% of each distribution. Shown with a 
diamond is the observed difference for the Milky Way, assuming a 
Ga lactic m etallicity of [Fe/H\ ~ —0.7 ±0.1 solar, as determined 
by llvezi~c t al. (20 0^, with a n a-enrichment of ~ +0.3 as sug- 
gested by . Chang et al l (Eooi) yielding Zuost ~ -0.4 ± 0.1. The 
stellar halo metallicity is taken to be Zjjji^ ~ —1.1 as discussed in 
§[4] The upper axis shows the central galaxy metallicities for each 
decade in host halo mass, as drawn from the fiducial Af^o^j — Af* 
relation and the case 1 mass-metallicity model. 



metal-poor intrahalo component from younger, brighter and 
more metal-rich stellar substructure currently undergoing 
disruption. In the left panel of Figure O we plot our case 
1 IHL metallicity as a function of host mass for only those 
subhalos whose tidal streams have dimmed to a B-band sur- 
face brightness lower than /is = 28 mag/arcsec^, and in the 
center panel we show the same trend for diffuse features dim- 
mer than ~ 27 mag/arcsec^. By contrast, the right panel 
excludes faint subhalo debris, selecting dynamically young 
streams that remain brighter than /ib ~ 27 mag/arcsec^, 
and thus becoming noisy at the galactic mass scale, since 
not every realization in this regime has had such an event. 



4 DISCUSSION 

One explanation for the origin of the mass-metallicity 
relation for star-forming galaxies arises from galactic 
winds that induce a significant amount of metal-enhanced 
mass l oss from galaxies with stellar mass below roughly 
j^qIo.s jTremonti et all |2004 ICallazzi etall |2005| , see also 
iKau ffmann e t al J regarding this threshold's importance 
for galaxy evolution in general), with the "blowout" process 
becoming more efBc i ent as the system's host mass decreases, 
although I Lee et al.1 (|2006h endorse a more leisurely model 
for metal outflows due to the small variance they obtain in 
the Mi, — Z relationship for dwarf irregular galaxies. Below 



the critical scale M^'' 10^°-^ Mq, the M* - Z correla- 
tion's slope grows steeper with lookback time, supporting 
a galactic chemical evolution model predicated on longer 
e- folding times for star formation in lower-mass galaxies 
(|Savaglio et al.ll2005l ). a hypothesis consistent with the em- 
pirical phenomenon of cosmic "downsizing" in which more 
massive galaxies form the majority of t heir stars at a higher 
reds hift than do less massive systems (|Heavens et al. I l2004l : 
Juneau et al]|2005l : Ferez-Gonzalez et allbOOTl i 

Numerical hydrodynamical simulations of hierarchical 
structure formation have also generated a characteristic stel- 
lar mass M^"^ ~ 3 X 10^" Mq which divides galaxy popula- 
tions into two groups. The first contains large systems that 
formed most of their stars more than ~ 10 Gyr ago, having 
been assembled via mergers of stellar-dominated subhalos 
with little gas left over to produce a significant starburst 
during coalescence. The second group is composed of com- 
paratively small systems that form stars more passively, al- 
lowing gas-rich mergers to produce enough star formation 
that the metallicity of the overall stell ar content is affected 
by the s ignificant fraction o f new stars (|de Rossi et al. I l2007l . 
see also lTissera et allbOOSl ). This stellar mass threshold re- 
mains constant to a redshift of z ~ 3, and M*'"** galaxies 
have roughly solar abundances that evolve very weakly over 
the same timeframe, these systems being typically embedded 
in halo s of mass Mhoat ^ 10 ^' ^'^ Mp , according to the formal- 
ism of Ivan den Bosch et al.l l|2007l ). Convolving this expec- 
tation with our prediction regarding the metallicity of dif- 
fuse intrahalo light prod uced by g alaxy disruption, we note 
that both simulations IStewart et al. 2003) and analytic 



models of galaxy formation ( Purcell et al] |2007| ) describe 
ACDM structure on all scales as being preferentially formed 
via mergers involving subhalos of mass M ~ G.lMhoat, in- 
dicating that both accreted substructure and the diffuse in- 
trahalo component should have roughly solar metallicities 
in halos of mass Mhoat ~ lO^'^'^M©. This context slightly 
favors our case 1 result, in which the mean metallicity of in- 
tragroup light on the relevant scale is very close to the solar 
oxygen abundance, although there is a scatter o ~ 0.15 dex 
at fixed host mass in both cases. 

Recent measurements of stellar halo metallicity in the 
Local Group se em to cl uster around the point Z~ —1.5 
(|Fergusonll2007l . see also iKalirai et al]l2006l for a review of 
the halo abundance of M31). However, the alpha-enrichment 
of the Gal actic halo is signi ficant (observed at -1-0.4 as re- 
viewed bv lMcWilliamlll997l '). which is consistent with an 
intrahalo stellar population forming in an early burst, be- 
fore t ype la supernova e co uld contribute m uch iron to the 
IGM (|Matteuccil |2003| . see IVenn et al. I l2004l for a review of 
chemical signatures in hierarchical formation models). An 
appropriate observational benchmark for the stellar halo 
oxygen abundance is therefore [0/H]/hl ~ —1.1 with re- 
spect to the solar metallicity, which falls along the lower 
68% contour in case 1 of our fiducial theoretical result. For 
the Galactic disk's iron abundance of \Fe/H\ ~ —0.7 ± 0.1 
with respect to the solar metallicity jlvezic et aL 20081). and 
an Q-enrichment of ~ -1-0.3 as suggested by Ichang et al] 
( 2002), we also note that the difference between this Zmw ~ 
—0.4 ±0.1 and the IHL metallicity is ~ 0.7 dex. We show 
the expected behavior of the quantity Zhoat — Zihl as a 
function of host mass in Figure |6] and note that the Milky 
Way's observed value slightly favors our case 1 model; fu- 
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Figure 7. Case 1 [0/-ff]iHL ^ function of IHL fraction for each realization at various fixed host masses. 



ture abundance surveys of distant galaxies and diffuse light 
may provide further empirical constraints supporting the 
prediction that intracluster stars have similar metallicities 
to that of the cluster's central galaxy, whereas intrahalo pop- 
ulations on smaller scales are significantly more metal-poor 
than their host galaxies. 

The evolution of a galaxy's metallicity arises from the 
equilibrating processes of the infall of gas and its subsequent 
enrichment in metals via star formation; these competing 
rates of dilution and enhancement are governed by the ef- 
ficiency of gas accretion and processing. Recent theoretical 
investigations have reproduced observed metallicities in the 
intergalactic medium at 2 > 2 and matched the — Z rela- 
tion's slope and normalization by introducing a paradigm in 
which galactic wind velocities are proportional to the stel- 
lar velocity dispersi on of the host system (t he "momentum- 
driven" scenario of iFinlator fc Pavel liooa ) and the mass- 
loading factor rj of the outflow scales with a galaxy's stellar 

— 1/3 

mass as 77 oc at low masses and approaches unity at 

the blowout point, where winds begin to escape their host 
halo and the M* — Z relation flattens to a roughly constant 
value. The steep rise in IHL abundance at and below the 
galactic scale occurs as intragalactic metal depletion pro- 
cesses dramatically fall to minimal efficiency with increasing 
host halo mass, thereby enriching satellite systems prior to 
disruption, and thus increasing the metallicity of the diffuse 
stellar component in higher-order structure as a similarly 
fast function of host mass. Given the rapidity of this trend, 
and the significant variance we obtain in the metallicities of 



stellar halos belonging to spiral-sized galaxies, we note that 
observational constraints on a possible correlation between 
IHL metallicity and host galaxy luminosity will likely reflect 
the stochastic and ongoing nature of the diffuse stellar com- 
ponent's formation, as well as the particular mass accretion 
history of a given galaxy. 

One such relationship is advocated bv iMouhcine et aP 

in which stellar halo abundances are derived from 
color measurements of the red giant branch in off-disk flelds 
and shown to increase with gala ctic lumin o sity a cross a wide 
range -2<Z<- 0.5, although iFergusonI (|2007l ') claims this 
correlation is non-existent when halo stars are kinematically 
selected. A more recent application of this method to an 
RGB-star survey of the Galaxy-analogue system NGC 891 
has shown a broad IHL abundance distribution centered 
around Z = —0.9, more metal-rich by a factor of three 
than the stellar halos of both large Local Group galaxies 
iMouhcine et al.ll2007l ). We do predict a strong rise in halo 
metallicity as host mass increases through the galactic scale 
(shown in Figure |4]), but the wide scatter at fixed mass sug- 
gests that the statistical accuracy of any such empirically- 
derived relation would be dramatically contaminated by 
variance in the individual mass assembly processes for each 
galaxy, as demonstrated by Figure [7] in which we obtain 
a power-law relationship between IHL abundance and the 
fraction of a system's stellar mass contributed by the in- 
trahalo component (n ~ 0.95 in [0/H]mL (x n log /ihl), 
spanning the widest range at galactic scales and thereby sug- 
gesting that stellar halo metallicity more properly correlates 
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with the luminosity of the stellar halo, and not necessarily 
the luminosity of the parent galaxy. 

This correlation grows weaker when we consider increas- 
ingly more massive host halos, as the IHL metallicity distri- 
bution tightens and the IHL fraction flattens. We predict an 
[O/IfJiHL value of ~ —0.5 for weak groups such as the Local 
G roup, which has a mass Mhoat ~ 5.3 x lO^'^M© according 
to lLi fc White ' (2008), although observational constraints in 
this regime remain elusive. On still larger scales, the Virgo 
cluster provides a more useful laboratory for the study of in- 
trahalo stellar metallicity, although systematic errors in the 
photometry of RGB-tip stars in these fields ca use great vari- 
ance in the abundance values yielded [e.g., iDurrell et al.l 
|2002| . in which the authors ob t ain - 0.8 < Z < -0.2). In 
a recent study, IWilliams et all (|2007h identify ~ 5300 in- 
tracluster stars in Virgo, determining that 70% — 80% of 
them are older than ~ 10 Gyr and populated over a wide 
range in abundance (— 2.5 < Z < 0.0), as we might expect 
from a diffuse population with multiple origins; the remain- 
ing 20% — 30% of the intrahalo population has a metal- 
licity distribution function peaking near solar values, these 
younger stars having been more recently orphaned by dis- 
ruption processes acting on their parent galaxy. Still more 
metal-rich is the intrahalo population of the cluster Abell 
3888, which contains one component having formed at red- 
shift z > 1 with high metallicity (1 < Z<2.5), as well as a 
more centralized d iffuse region at roughly solar abundance 
iKrick et al.ll2006l l. The variation in our predicted distribu- 
tion is of order o" ~ 0.1 dex, although the location of the 
plateau varies significantly with surface brightness cuts in- 
tended to separate older intrahalo populations from dynami- 
cally young IHL plumes in galaxy clusters. Typical measure- 
ments of intracluster light have analyzed fields without obvi- 
ous luminous substructure, corresponding to the spheroidal 
halo of field stars in spiral galaxies. It appears that we must 
be careful about claims as to how much stellar mass actu- 
ally lies in long-disrupted satellite remnants, as opposed to 
relatively bright diffuse tidal streams. 

However, it is worth noting that stellar halo observa- 
tions of galaxies outside the Local Group cannot distinguish 
stream-like substructure from background field stars, just as 
intragroup and intracluster studies cannot. Tidal features 
such as M31's giant southern stream are rich in stellar mass 
(the progenitor may have been as bright as a billion solar 
luminosities according to iFont et aL Il2006l) as well as met- 
als {Z ~ —0.5 according to Guhathakurta et al.ll2005l ). so 
we can reconcile our predictions with current galactic mea- 
surements given the caveat that typical halo metallicities 
reflect the underlying stellar component, without the inclu- 
sion of bright outer-halo substructure resulting from mas- 
sive accretion events, a distinction made more clear by a 
recent semi-analytic treatment of baryon dynamic s in nu- 
merical models of L^-galax y halos l|Font et al.ll2007l . see also 
iBullock fc Johnstonll2005l for details of the simulation). We 
also predict a clear correlation between the oxygen abun- 
dance and the surface brightness of tidal streams in galaxy 
clusters, a prediction that will be testable as future deep 
surveys come online and metallicity distribution functions 
become more well-studied in diffuse intrahalo stellar popu- 
lations at and above the galaxy scale. 
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